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c-| ! It is obtained that dark energy emerges from the higher-derivative 

^ | gravity with non-linear terms proportional to R 2 and R 3 with R being 



the Ricci scalar curvature. Interestingly, it is found that the universe 
begins with acceleration, which continues for a short period. Later on, 



it decelerates driven by radiation and subsequently by matter. Two 
components of dark energy are found here. In the early universe, 
dark energy mimics quintessence giving high value of initial density 
~ 10 75 GeV 4 at Planck scale. But, in the late universe, dark energy be- 
haves like phantom giving current value of its density and causing late 
acceleration for some time. Thereafter, deceleration driven by mat- 
ter resumes. Thus, two components of dark energy are obtained here 
(i) quintessence-like in the early universe and (ii) phantom-like in the 
late universe. Moreover, it is interesting to see acceleration being a 

transient phenomonon in early as well as late universe. 
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Modifying earlier results [1], observations of 16 Type SNe la show a 
conclusive evidence for deceleration preceding acceleration in the late 
universe [2] . Cosmic acceleration is fuelled by a source of energy violat- 
ing strong (SEC) or weak energy condition (WEC), called dark energy 
(DE). Instead of the simplest A- model (A being the cosmological con- 
stant) having fine-tuning problem, dynamical models, given by differ- 
ent scalars (quintessence, tachyon and phantom) were envisaged in the 
past few years to explain late cosmic acceleration and fall of DE den- 
sity (DED) from a very high to current low value [3]. Due to negative 
pressure and being better fit with WMAP, generalized Chaplygin gas 
(GCG) have been a good candidate for DE [4]. In [5], another attempt 
is made adding string gravity corrections with dilaton or modulous 
fields to gravitational action. Recently, some gravitational models are 
proposed, where non-linear term of the Ricci scalar R is taken as DE 
and its various cosmological consequences are discussed [6]. 

In these models, non-gravitational or gravitational source of DE is 
taken a priori. But, it is more interesting to have model, where DE 
emerges spontaneously from the gravitational sector. 

In what follows,, non-linear term of curvature is not taken as DE 
lagrangian a priori, which is unlike the approach in [6] . Rather, in 
the modified Friedmann equations (MFE), obtained here, DE density 
terms emerge due to presence of linear as well as non-linear terms of 
curvature scalar in the action. This is the approach adapted in [7]. 

In modern cosmology, universe is supposed to have acceleration in the 
early universe followed by deceleration driven by radiation and subse- 
quently by matter as well as late acceleration driven by DE, preferably 
phantom-like. Here, a cosmological model is proposed yielding this 
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cosmic scenario. In addition to DE terms (mentioned above), here, it 
is found that dark matter (DM) also emerges from the gravitational 
sector. For the radiation term another scalar <p is introduced. It is 
found that this term emerges as a combined effect of and the scalar 
curvature R. 

Here, it is found that DE is induced by curvature [7] causing early 
and late cosmic acceleration, which are transient phenomena. Inter- 
estingly, DE, obtained here, has two components (i) quintessence-like 
giving high value of DED in the early universe and (ii)phantom-like in 
the late universe giving current value of DED. This model is interesting, 
in the sense, that it explains different phases of the cosmic dynamics 
i.e. (i) inflation in the early universe, (ii) deceleration after exit from 
inflationary phase, (hi) transition from deceleration to acceleration in 
the late universe and (iv) another transition from acceleration to de- 
celeration in the future universe. Natural units (h — c — 1) are used 
here with GeV as a fundamebtal unit and lGeV -1 = 6.54 x 10 _25 sec. 

As mentioned above, modifications of gravitational action with higher 
powers of scalar curvature R is not new. So it is natural to probe cosmo- 
logical consequences on taking R 3 in the gravitational action. Earlier, 
Amendola [8] also had taken gravitational action having lagrangian 
R + R 3 . Here, investigations begin from the Planck scale. 

The action is taken as 



where G = Mp 2 (M P = 10 19 GeV is the Planck mass), a is a dimen- 
sionless coupling constant, (3 is a constant having dimension (mass) -2 
(as R has mass dimension 2)and V M stands for covariant derivative. 



S = 



I 




r R 



fai? 2 -^ 3 + {-^V M 0(x)V i ,0(x)-\/(0( a ;))} , 

(1) 



L16ttG 
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Moreover, <f)(x) is a scalar field with V((f>) as potential giving radiation 
term in the modified Friedmann equation obtained here. 
Action (1) yields field equations 

Y^q( R ^ - \a^R) + a{2V^V u R - Ig^UR - ^ U R 2 + 2RR^) 
-(3{3V^ V R 2 - 3 9txv UR 2 - X -g^R? + 3R 2 R^) 

+\{v il <b{x)V u <b{x) - g^(U°<P{x)V a <P{x) - = (2a) 

using the condition dS/dg^ = 0. Moreover, invariance of S under 
transformation <p — > <p + 5<p yields 

□0 + V'(<t>) = 0, (26) 

where V'{<pi) = dV(4>)/d(f> and the operator 

1 d ( , „„ d 



□ = 



-g dx^ 

Taking trace of (2a), it is obtained that 
R 



(S=»<r w ). (2c) 



167rG 



Q(a - 3(3R)DR + 18/3 ViJV^i? - (3R 



3 



+^[-V>W V CT 0(x) + 4\/(0(x))] = (3) 

In (3), (a — 3(3 R) emerges as a coefficient of DR due to presence of 
terms aR 2 and (3R 3 in the action (1). If a = 0, effect of i? 2 vanishes and 
effect of i? 3 is switched off for (3 = 0. So, an effective scalar curvature 
i? is defined as 

1 R = a - 3(3R, (4a) 

where 7 is a constant having dimension (mass) -2 being used for di- 
mensional correction. For t > t P (t P being the Planck time), scalar 
curvature at Planck scale Rp > R. Here a and (3 are taken so that 
(a — 3(3R P ) > 0, giving 

(a - 3(3R) >(a- 3(3R P ) > 0. (46) 
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Connecting (3) and (4a), it is obtained that 

.,2, 1 



R 



67 



1 a' 
+ 



167rG 3/5 



67 



1 a' 
+ 



167rG 9/3. 



7.R 



+ ^[7i? - 3a] + il-VXxjV^ix) + 4\/(0(*))] (5) 



Energy-momentum tensor components for are 

7J" = V^V^ - <T (±V>V^ - V(0)) 



yielding energy density for as 



and pressure for as 



p = ^V CT 0V CT + \/(0) 



I V CT 0v CT 0-\/(0). 



(6a) 



(66) 



(6c) 



Experimental evidences [9] support spatially homogeneous flat model 
of the universe 

dS 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 ] (7) 

with a(t) being the scale-factor. 

In this space-time, (2b), (5), (6b) and (6c) are obtained as 



a ■ 



+ 3-0 + V'(<f>) = o, 
a 



(80.) 



a R 67 
R 



1 a 2 
+ 



16ttG 3(3 
54/5 47^ 



1 

67 



1 a 2 
+ 



167rG 9/3 J 7j R 



2 + 4K(0) , 



and 



= \tf + v(4) 



p<t> = \<i> 2 - v(4) 



(86) 
(8c) 

(8d) 



as 0(x, i) = 0(t) due to spatial homogeneity. 
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For a(t), being the power-law function of cosmic time, R ~ a~ n . 
For example, R ~ a~ 4 for radiation model and R ~ a~ 3 for matter- 
dominated model. So, there is no harm in taking 

R=—, (9) 

where n > is a real number and A is a constant with mass dimension 
2. 

(8a) integrates to 



d 
dt 



1 12 



2 + U(0) +3-0 2 = O. (10a) 



L2 

For a cosmic fluid, pressure p is proportional to its energy density p. 
So, for fluid of scalar 0, the equation of state is taken as 

V<t> = k P<t» 

where k is a proportionality constant. Also, from (8c) and (8d), we 
obtain 

1 (l-k\ 
+ 

In what follows, we find very intersting results, on taking k suitably, 

as 

k=^. (10c) 
Using these values of k in (10a), we obtain 



Using (10c) in (10a) and integrating, we obtain 



Here B is an integration constant. 
From (lOd) and (11) 
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+ w(<j>) = 



3n 



B 



(2 + n) a( 2 "+ 4 ) ' 
Using (9) and (12) in (8b) ,it is obtained that 



(12) 



d /a\ ,„ n . /a\ 2 Ca n 



Da r 



+ 



7 A r 1 



54nf3 



3a 



3/3B 



where 



and 



C 



D 



4(2 + n) 7 3 A 2 a 4 ' 
1 



1 

67 

1 

67 



+ 



a 



16nG 3(31 



1 a^ 
IQ71G + 9(3 



(13a) can be re- written as 

Ca n+l 



a + (2 - 2n) 



d 2 



r Da"" 

I 1 " c. 




- 1 


3a " 


+ 54n/3 


.a" 


7A 



3/35 



which integrates to 



4(2 + n) 7 3 A 2 a 3 : 
2C 



i 6 " 4 " + nA 



Da r 



>yA 



(6-3ra) C(6-2ra)] 
1 3a 



27n/5 L (6 - 5n)a" 7^(6 - An) 
3(3B 



(13a) 

(136) 
(13c) 



(14) 



15a) 



2(2 + n)(2-4n) 7 3 A 2 a 4 ' 
where E is an integration constant having dimension (mass) 2 . 

(15a) is the modified Friedmann equation (MFE) giving cosmic dy- 
namics. In this equation, the term proportional to a -4 is the combined 
effect of curvature and scalar field 0. Other terms, on r.h.s. of this 
equation, are imprints of curvature. The first term of these, propor- 
tional to a~( 6 ~ 4 ") emerge spontaneously. It is interesting to see that 
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this term (the first term on r.h.s. of (15a)) corresponds to matter den- 
sity if n — 3/4 i.e. for this value of n, it reduces to Ea~ 3 . As it is 
spontaneously created from the gravitational sector, we recognize it as 
as dark matter density. 

For n = 3/4, (15a) looks like 

5,0a 3 / 4 1 



E 32Ca 3 / 4 r 



a° 



I67A 



45A 
1 



729/? 



,3/4 



1 - 

a 



6C J 



Q/3B 



ll 7 3 A 2 a 4 



(15b) 



On r.h.s. of (15b), terms proportional to a 3//4 (a~ 3 / 4 ) are recognized 
as dark energy density being created by curvature. Moreover, the term 
proportional to a -4 corresponds to radiation as usual. 

(15b) is re- written as 



32Ca 3 / 4 r 5Da 3 / 4 



45A 



1 - 



6C 



I67A 
729(3 



,3/4 



a 

>yA 



57T 



G 



with 



and 



Prd 



9(3B 



3E 



(Prd + Pdm) 

(15c) 
(15d) 



Pdm = 8^" (15e) 
Using current values of p r d and pd m as 5 x 10~ 5 p cr and 0.23p cr re- 
spectively [10], (15d) and (15e) yield 



Pad 



0.23p c 



Oo\3 

a ) 



and 



Prd 



5 x 10 >,.,.' 



where a = a(f ), 3E/SnG = 0.23p cr a^, 9(3B/AA-f 3 A 2 
and 

_3Hl 



(16a) 

(166) 
5 x 10~ 5 p cr a 4 
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with Hq being the current Hubble's rate. 
(16a) and (16b) look like 



and 



0-23p cr 

Pdm = 5— (16c) 

a 6 



5 x 10 5 p cr 
Prd = 5 (16d) 

(X 



on normalizing ao as 

a = 1- (16e) 
(16c) and (16d) yield that p rd < p m at red-shift 

Z=--l< 5399. 
a 

WMAP [10] gives decoupling of matter from radiation at red-shift 

Z d = — - 1 = 1089. (17) 
It shows that dark matter dominates over radiation when a > ad- 
Case 1. When a < 2.17 x IP' 4 

This is the case of early universe, where (15c) reduces to 

er^[{^--i+ ftd ] (18) 

as term ~ a 3 / 4 dominates over term ~ a 3 / 4 . In (18), 

_ 16 7 A 3/4 



and 



A qu = « / 196 n 

243ttG/3 1 J 

being the cosmic tension in the early universe, which is caused by scalar 
curvature. 

The conservation equation for DE is 

p de + 3-(l + w de )p dc = 0, (20a) 

CI 
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where equation of state parameter (EOSP) Wd c = Pde/pde with p^ e 
being DE fluid pressure . (19a) and (20a) yield 

< = ~ (206) 

It shows that DE component, dominating the early universe, mimics 
quintessence. 

At the Planck scale, pde = (3/8irG)Mp, giving 



qu 
Pdc 



8ttG 



M 



/a P \3/ 4 



(21) 



with ap = a(tp) and tp = M p l being the Planck time. 
Using (16d), (19a) and (19a) in (18), it is obtained that 



I67A 



ay 

a) " 729/3a 3 / 4 



a 



1- — -a 



3/4 



+ 



5 x 10~ 5 H$ 



(22a) 



(226) 



This equation shows that at 

a = a qe =(— )j , 

terms within the bracket on r.h.s. vanishe and, for a > a qe , expansion 
of the universe is driven by radiation. 

In the very early universe, first term on r.h.s of (22a) dominates over 
the third term, so (22a) is obtained as 



ay 



I67A 



a J 729/3a 3 / 4 



1 - ^-a 3 / 4 



(22c) 



which integrates to 



aa 



3/4 



a 3 / 8 = a 3 / 8 (^^) 1 sin [V27rGAi u (t - t P ) + sin' 1 ^ y — 

As, in the early universe, (t — t P ) is very small and \J aa 3 / 4 ~J%A is also 
small, so this solution is apprximated as 

8/3 



a(t) ~ a P ^(27vG-fA/a)\i u (t -t P ) + l 
using sinx ~ x for small x. 



(22d) 
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Also, at Planck scale, p^ (p) > 2A qu . So, from eqs.(19a), (19b), (22b) 
and inequality (4b) 



3/3R P < a < 



87A 



3a 



3/4" 



Connecting (22b) and (22d) 



8/3 



(22e) 



a(t) ~ a P ^/27rGA ( i u (a q e/ap) 3 /4(t - t P ) + 1 

(22e) gives power-law inflation (acceleration) in the early universe. 
This scenario continues till p*^ = 2A qu upto the time 



t P + 



as 



3/4 



(23a) 
(236) 



By the time t qe the scale factor grows to a qe . If the early universe 



inflates sufficiently during this period and goes over to 65 e-foldings, 



^qe ^ 10 28 
dp 



(23c) 



and 



tqc — tp 



2t, 



,y|^) 3/8 = 5.I6 x 10 10 tp = 3.37 x lCT 27 sec. (23d) 



3 V ap J 

Connecting (23b) and (23c), it is obtained that 



A qu 



-M 



2^ 3/4 = 5. 97 xl0 53 GeV 4 



167rG 

For t d >t>t qe , using = 2A qu (22a) is obtained as 

2 „ /1\ 4 



(23e) 



which integrates to 

a(t) = aJl + 0.UU2H o (—)\t-t qe ) 
l \ a qe / 

giving deceleration as a < 0. 



1/2 



(24) 
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Case 2. When a > 2.17 x 10~ 4 

This is the case of late universe, where (15a) reduces to 

as term ~ a 3 / 4 dominates over term ~ a -3 / 4 . Here 



and 



ir 2 

Equation, like (25) with p 2 term, also arise in brane-gravity inspired 
Friedmann equation with A called as brane-tension [11]. Here, this 
term arises from the 4-dimensional higher-derivative gravity without 
using any prescription of the brane-world. A ph is called cosmic tension 
in the late universe. 

Connecting (26) and (20a) 

w d P c h = -5/4, (28) 

which shows that DE component dominant in the late universe mimics 
phantom. 

Current value of pde is 0.73p cr [10], so (26) reduces to 

pt = 0.73 Pcr a 3 / 4 (29) 

with AC/lhuAG = 0.73p cr . For a < a = l,pg > (pS) 2 - So > for 
Pdm > Pdc, (25) integrates to 

a(t) = a d [l + 25887.8F (t - t d )] 2/3 (30) 

for (l/a d ), given by (17) and p dm given by (16a). (30) gives a < 
showing deceleration. 
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(16a) and (26) show transition from p m > p^ 1 to p m < p^ at red-shift 

1 /73\ 4 / 15 
Z. = _-l=(_) -1 = 0.36 (31) 

with a* = a(t*) and t* being the transition time. Thus, red-shift, corresponding 

to the transition time £*, is obtained within the range 0.33 < < 0.59 

given by 16 Type supernova observations [2]. 

Thus, for t > £*, using (26) and (28), (25) reduces to 

/a\2 _ SnG h f _ pf e \ 
\a) ~ 3 Pdc \ 2AP h J' 

which integrates to 



a(t) 



0.73 Pcr r / 0.73 Pc: 
2A^ + \V L26 "^ 

2-i -4/3 

^ H<>\/{).7-M1 - f , ) " 



-^ov/oTrs^-t*)} 2 ] 



(32) 



using (29) and (31). 

(32) shows acceleration till p^ becomes equal to 2A ph as it grows 
with expansion. It happens till a(t) increases to a pc satisfying 

^ - Sr„ < 33 > 

Thus, accelerated growth of a(t), given by (32), stops at time 



t pe = U + ^== \l 1.26 - (34) 

P 3 J ff v / a73V 2A? h 

and (25) reduces to 

/d\2 8ttG 



yielding 



a(f) = a pe [l + 0.72i/ a- 3 / 2 (t - t pc )] 2 / 3 (36) 



Thus, late acceleration is obtained during the period (t pe — t*). For 
t > t pe , DE is not able to drive cosmic dynamics as p^ > 2A ph leads to 
DED component negative in (25) and a(t) complex. So deceleration, 
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driven by dark matter, will resume. Interestingly, solution (32) does 
not fall into finite time future singularity [12]. Here, late transient 
cosmic acceleration is obtained using the modified gravity. In [13], this 
phenomenon is obtained using RS2 brane-model. . 

Thus, it is found that DE, emerging due to curvature, causes tran- 
sient acceleration in the early and late universe . It is transient as 
cosmic tension puts a brake on accelerated expansion. Interestingly, 
in early and late universe, cosmic tension also arises from the grav- 
itational sector. Moreover, it is interesting to see that quintessence 
dark energy emerges in the very early universe, where scalar curva- 
ture is high and phantom dark energy emerges in the late universe, 
where scalar curvature is low. Thus, according to the present model, 
low scalar curvature produces phantom dark energy and high curvature 
produces quintessence dark energy. 
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